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Abstract. Polygons are described as almost-convex if their perimeter differs from 
' the perimeter of their minimum bounding rectangle by twice their 'concavity index', 

• to. Such polygons are called m-convex polygons and are characterised by having 

fH | up to m indentations in their perimeter. We first describe how we conjectured the 

(isotropic) generating function for the case m = 2 using a numerical procedure based 
on series expansions. We then proceed to prove this result for the more general case 
of the full anisotropic generating function, in which steps in the x and y direction are 
distinguished. In so doing, we develop tools that would allow for the case to > 2 to be 
^ , studied. 
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1. Introduction 

The enumeration of self-avoiding polygons (SAPs) is a classical problem in statistical 
mechanics and combinatorics. Exact results have thus far largely required the restriction 
of SAPs to sub-classes that are in some way convex. In two dimensions, convexity means 
that the perimeter is equal in length to the length of the minimum bounding rectangle 
(MBR). Column- convexity means that any vertical cross-section may only intersect 
the polygon twice, such that all columns are connected. Examples of convex and 
column-convex polygons can be seen in figure 1. Convex polygons on two-dimensional 
lattices have been studied extensively by Lin [1, 2, 3] and Bousquet-Melou [4] and 
many exact results are known including the full area-perimeter generating function. In 
1997, Bousquet-Melou and Guttmann [5] gave exact results for convex polygons in three 
dimensions and a method for their enumeration in arbitrary dimensions [5]. 



(a) Pyramid. (b) Unimodal. (c) Convex. (d) Column-convex. 

Figure 1. Examples of polygons that are in some way convex. 

Enting et al [6] described polygons as almost-convex if their perimeter differs from 
the perimeter of their minimum bounding rectangle by twice their 'concavity index', 
m. Such polygons are called m-convex polygons and are characterised by having up to 
m indentations in their perimeter. Examples of 1-convex and 2-convex polygons can 
be found in figure 2. Enting et al derived the asymptotic behaviour of the number 
of m-convex polygons according to their perimeter, n for m = o( v m~). The results 
were confirmed for the case m = (i.e. convex polygons) by the known perimeter 
generating function. Subsequently, Lin [7] derived the exact generating function for 
1-convex polygons, using a 'divide and conquer' technique introduced to the problem 
of convex animals (the interior of a convex SAP) by Klarner and Rivest [8]. His result 
provided support for a conjecture in [6], giving the next term in the asymptotic expansion 
for the number of polygons with perimeter n and concavity index m. 

This is the second in a series of papers that look at families of m-convex polygons. 
In the first [9] we outlined the 50-year history of polygon enumeration on the square 
lattice before re-deriving the generating functions for 1-convex polygons in an effort 
to generalise the methodology and extend the results to osculatingf and neighbour- 
avoiding^: polygons. 

f Osculating polygons are those that may touch themselves, but not cross. 

\ Neighbour-avoiding polygons are those that may not occupy a neighbouring lattice vertex without 
being connected by an edge. 
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(a) A 1-convex polygon. (b) A 2-convex polygon. 

Figure 2. Almost-convex polygons with the minimum bounding rectangle marked. 

Polygon models have long been used to model vesicles, with self-avoiding polygons 
being the canonical model [10] . Associating a fugacity with the area of the polygon, a 
phase transition occurs, so that at sufficiently high fugacity, the polygons change from 
the self-avoiding class to the convex class. The introduction of m-convex polygons 
permits the exploration of this phase transition in greater detail. Firstly, we find 
the (unsurprising) result that m-convex polygons, for m finite, have the same fractal 
dimension (and hence are in the same universality class) as convex polygons, that is to 
say, the fractal dimension remains unchanged at 2 (as compared to the value 4/3 for self- 
avoiding polygons). Secondly, the m-convex model permits one to associate a fugacity 
with the concavity index m, and this would be the polygon analogue of the stiffness 
in self-avoiding walk models of polymer stretching. It is however not our purpose to 
study this aspect of the problem here. Ideally, one would like to predict how the form 
of the generating functions changes as the concavity index grows. This would help us 
understand what happens in the scaling limit as the concavity index grows in proportion 
to the perimeter. 

In 2000, in unpublished work, we conjectured, on the basis of long series expansions, 
the (isotropic) generating function for the case m = 2. In this paper, we describe the 
numerical procedure that led to this conjecture, and proceed to prove it. Indeed, we do 
so for the full anisotropic generating function, in which steps in the x and y directions 
are distinguished. In so doing, we develop tools that allow for the case m > 2 to be 
studied, though we do not do so. 

In section 3 we describe the series expansions that allowed us to conjecture the 
exact result, and also, not incidentally, to provide checks on our rigorous results in the 
process of proving our conjecture. In the following section we introduce the methodology 
used to derive our results, followed by examples of its application. We enumerate all 
the separate building blocks required in the factorisation of 2-convex polygons. The 
intermediate results, as well as much of the detail, are omitted for reasons of conciseness 
and clarity. These may be found together with the equivalent staircase and unimodal 
results in [11] and [12], where the presented results first appeared. 
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Figure 3. A convex polygon with marked vertices where the polygon meets the 
minimum bounding rectangle 



2. Definitions and notation 

In this section we briefly summarise some definitions and notation used in the remainder 
of the paper. 

Directed Walks (DWs): Square lattice walks that take either positive or negative steps 
in each of the horizontal and vertical directions, but not both. For example, a walk 
that only steps up and to the left is directed. 

Generating functions: If the number of polygons (in a given class) with perimeter n 
is p n then the associated isotropic generating function is F(x) = J2 n Pn xH - F° r 
polygons on the square lattice the perimeter is always even and we shall therefore 
study the half-perimeter generating functions G(x) = ^2 n P2 n x n - In more generality 
we distinguish between the number of steps in the x and y directions and study 
the full anisotropic generation functions, G(x,y) = J2 m ,nP^m,2nX m y n , where p 2m ,2n 
is the number of polygons with 2m horizontal steps and 2n vertical steps. 

The half-perimeter operator: We denote by E the operator that converts the perimeter 
generating function to the half-perimeter generating function. E x (resp. E y ) 
converts only the direction counted by x (resp. y). (For a full definition, see 
[5].) We therefore have 

E[f(x,y)}=E x [E y [f(x,y)]] (1) 

where 

E x [f(x, y)] = (/(Vi, y) + f(-VZ, y))/2. (2) 

If x or y (or a function of x or y) is asterisked, then the operator only takes the 
half-perimeter of the non-asterisked generating function. That is, 

E x [f(x,x*)]=E x [f(x,y)]\ y ^ x . (3) 

For example, E x [l/(1 — x — x*)] — (1 — x)/(l — 3x + x 2 ) and, more generally, 

E x [f(x)g(x*)]=g(x)E x [f(x)}, (4) 
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which extends naturally to the multivariate case. This allows the conversion from 
the perimeter generating function to the half-perimeter generating function for 
different factors of a polygon separately. 

Minimum bounding rectangle (MBR): This is the smallest rectangle which encloses the 
polygon. The polygons in figure 2 have a 6 x 6 MBR while the polygon in figure 3 
has a 12 x 8 MBR. 

Sides: A convex polygon can be factored, as per figure 3, into four overlapping DWs: 
from li to £2, from t\ to r±, and so on. These are the maximal directed factors of 
the polygon, and we refer to them as the sides of the polygon. 

Arcs: An arc of a convex polygon is made up of a pair of adjacent sides. That is, it 
a maximal partially-directed factor of the polygon. Referring to figure 3, the top 
arc is the path from l\ to r±, passing through ^2,^1, ^2 and r 2 , and similarly for the 
bottom, left, and right arcs. 

Indents: An indent occurs in a side when the DW takes a step in the 'wrong' direction. 
We refer to such walks as 'almost-directed'. As an example, consider the part of 
the perimeter on the top-right side, which, when the polygon is traversed anti- 
clockwise, only takes steps up and to the left. An indent would occur if this walk 
were to take steps down (or to the right) and then resume taking steps up and 
to the left. We shall distinguish between indents in the vertical and horizontal 
directions. In figure 2a the polygon has a single vertical indent while the polygon 
in figure. 2b has both a vertical and horizontal indent. A walk taking k' steps in the 
'wrong' direction followed by k" steps in the 'correct' direction produces an indent 
k = mm(k', k") deep. Note that a vertical indent may contain a horizontal indent 
and vice versa. 

3. Exact solution from enumerations 

Several years ago two of us (IJ and AJG) found the exact generating function, C 2 (x) for 
2-convex polygons numerically from exact enumerations for 2-convex polygons and some 
simpler sub-classes. We found the solution by counting the number of 2-convex polygons 
by using a program designed to enumerate SAPs [13]. This program counts the number 
of SAPs by perimeter having a given MBR. From this data it is trivial to extract the 2- 
convex data. However, with the computational resources available at the time, we could 
not get a series long enough to find the generating function directly (we counted 2-convex 
polygons up to perimeter 110 yielding 48 non-zero terms). We therefore enumerated 
three sub-classes, namely 2-convex polygons with one vertical indent 2 deep on the top 
arc, 2-convex polygons with two vertical indents 1 deep on the top arc, and 2-convex 
polygons with vertical indents 1 deep on the top and bottom arcs. The programs used 
in these enumerations were simple generalisations of the one used by Guttmann and 
Enting [14] to count convex polygons. We then used the series for these sub-classes to 
find the respective generating functions F 2 (x), F s (x) and F Q (x). Finally, we looked at 



Exact generating function for 2-convex polygons. 



6 



the 'remainder' of the full 2-convex case F R (x) = C 2 (x) — AF 2 (x) — 4F s (a;) — 2F (x), 
which counts cases of 2-convex polygons with vertical and horizontal indents 1 deep, 
and managed to find the generating function (and hence the full 2-convex generating 
function). 

We found the solution to the special cases via a judicious guess for the form of 
the generating functions. From Lin [7], the isotropic generating function for 1-convex 
polygons is: 



From this it is reasonable to expect that the generating function for 2-convex 
polygons C 2 (x) and the special cases F 2 (x) and so on are of a similar form, [A(x) + 
B{x)yJ\ — Ax]/D(x), where A(x), B(x) and D(x) are polynomials. In particular we 
expect the denominator D(x) to be similar to the one in the above expression, but 
with larger exponents and possibly involving further simple factors. If we can find 
D(x) then A(x) and B(x) can be found simply from a formal series expansion using 
polynomials with unknown coefficients. By equating the terms in this formal series 
with the known series for say F 2 (x) we find a set of linear equations for the unknown 
polynomial coefficients. 

Here we give some further details of how we found the generating function F 2 (x). 
We calculated the number of convex polygons with an indent 2 deep on the top arc to 
perimeter 216. This gives us the first 100 non-zero terms in the half-perimeter generating 
function. Our first task is to determine the denominator D(x). We did this by analysing 
the singularities of F 2 (x) using differential approximants. Our analysis showed that the 
series has singularities at x — 1/4 with exponents —3 and —5/2 (this confirms that at the 
dominant singularity we have a square-root correction term), at x — 0.381966. . . (the 
first root of 1 — 3x + x 2 ) with exponents —3 and — 1, at x — 1 with exponents —4.99(3) 
and —3.0(5), and at x — 2.618 . . . (the second root of 1— 3x+x 2 ) with exponent —3. The 
conclusion is that in this case D(x) = (1 — x) 5 (l — 3x + x 2 ) 3 (l — 4a;) 3 . By inserting this 
into the general form and equating terms in the formal expansion with those of F 2 (x) 
we found a solution with polynomials A(x) and B{x) of degree 18 and 17, respectively, 
with B(x) containing the factors (1 — x) 2 {\ — 3x + x 2 ) 2 , the latter of these factors was 
indicated by the exponents found at x — 0.381966 . . .. The polynomials are: 



Ci(x) 



x 3 (-4 + 56a; - 300a; 2 + 773a; 3 - 973a; 4 + 535a; 5 - 90a; 6 + 24a; 7 ) 
(l-a;)(l-3a; + a; 2 )(l-4a;) 3 
4a; 3 (1 - 9a; + 25a; 2 - 23a; 3 + 3a; 4 ) 

+ (l-a;)(l-4a;) 5/2 ' 



(5) 



A(x) 



B(x) 



-8a; 2 + 208a; 3 - 2428a; 4 + 16856a; 5 - 77742a; 6 + 252114a; 7 
-593563a; 8 + 1032521a; 9 - 1336471a; 10 + 1284072a; 11 - 904540a; 12 
+456064x 13 - 158327a; 14 + 36093a; 15 - 4955a; 16 + 126a; 17 + 88a; 18 . 

(1 - a;) 2 (l - 3a; + a; 2 ) 2 (8a; 2 - 128a; 3 + 844a; 4 - 2992a; 5 
+6262a; 6 - 8014a; 7 + 6188a; 8 - 2602a; 9 + 470a; 10 - 12a; 11 ). 



(6) 
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Similarly we found the generating functions for the other three special cases Fs(x), 
Fq(x) and F^(x). The only additional point worth noting is that the denominator in 
the case of -Fr(x) contains the extra factor (l — 2x). Collating these results we find that: 

s> / \ M x ) B 2 (x) 

2{ ' (l- X y(l-2x)(l-3x + x 2 ) 3 (l-4x) 4 (l-a;) 3 (l-4a;) 7 / 2 ' K ) 

where 

A 2 (x) = -24a; 2 + 864a; 3 - 14368a; 4 + 146672a; 5 - 1030216a; 6 + 5289512a; 7 
-20587766a; 8 + 62176564a; 9 - 147946110a; 10 + 280112802a; 11 
-424512212a; 12 + 516373058a; 13 - 504068274a; 14 + 393649476a; 15 
-244279626a; 16 + 119050550a; 17 - 44773540a; 18 + 12722814a; 19 
-2660520a; 20 + 378184a; 21 - 22560a; 22 - 3200a; 23 + 512a; 24 . 

B 2 (x) = -24a; 2 + 456a; 3 - 3592a; 4 + 15264a; 5 - 38200a; 6 + 57792a; 7 

-52832a; 8 + 28872a; 9 - 8968a; 10 + 1248a; 11 + 128a; 12 . (8) 

In the next sections we show how to prove this result for the general anisotropic 

case. 



4. Enumeration techniques 

4-1. Convex polygon basics 

Following from Section 2, one may describe convex polygons as a series of four non- 
intersecting DWs that make up the four sides of the polygon. (To see this, we refer to 
figure 3.) If one of these sides has no steps in the interior of the MBR, this means that 
the polygon touches one of the corners of the MBR. This sub-class of convex polygons 
is referred to as directed-convex§ or unimodal||. We could therefore define this class 
of polygons as having only three sides. Similarly, staircase polygons are formed of two 
DWs that start and end at diagonally opposite corners of the MBR. Finally, pyramids 
and stack polygons also have two sides formed by DWs, but these are adjacent sides, 
with a straight base or side edge. (Stack polygons are simply pyramids on their sides.) 
Such classes of convex polygons were depicted in figure 1. 

For notational convenience, let us now define some well-known generating functions, 
where x (resp. y) counts the horizontal (resp. vertical) steps. The generating function 
for pairs of intersecting DWs that begin and end at the same points (referred to as 
staircase festoons) we denote as 




§ This name comes from the definition which says that all cells in the interior of the polygon can be 
connected with the corner cell by a directed walk on the dual graph. 

|| This name comes from the fact that there is only one mode in each direction, when we take the 
projection of the walk in that direction. 
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where A = 1 — 2x — 2y — 2xy + x 2 + y 2 . The staircase polygon generating function we 
denote as 

S = (1 - x - y - VA)/2. (10) 

We note that the unimodal generating function is simply xyZ. Now, by defining 

u = x + S and v = y + S, (11) 

we can re-express all our almost-convex polygon generating functions as expressions 
with terms that are simply the quotient of polynomial functions of u and v. This is 
achieved via the transformation of variables 

x — u(l — v) and y — v(l — u). (12) 

For example, we have 

A = (l-u-v) 2 , Z = \/{\-u-v) and S = uv. (13) 

4-2. Joining polygon factors 

The Temperley method is central to the enumeration of partially convex polygons. 
The so-called 'functional- Temperley' method allowed Bousquet-Melou [4] to enumerate 
classes of column-convex polygons. It differs from the Temperley method in that 
it can be used to concatenate several large enumerable parts of the polygon, rather 
than individual columns. In particular the concatenated building blocks may be 
different types of polygons. A variation [15] allowed for the enumeration of certain 
classes of animals, represented as heaps of dimers. Rechnitzer [16] identified these 
methods as equivalent, the superiority of one over the other lying in its ease of use and 
appropriateness to the recurrence relation underlying the problem. 

We call distinct enumerable parts of the polygons factors, due to our ability to 
factorise the polygons into such parts by separating them at unique factorisation points. 
In this paper, we use factorisation lines, which we define by extending the interior edge 
of the indentations of almost-convex polygons into lines that bisect the lattice (see 
figure 8.) 

4-2.1. The Hadamard product The approach that has been used repeatedly by Lin 
[1, 2, 3] in the enumeration of convex polygons is to build them up vertically, block by 
block. The functional- Temperley method can therefore be used in this case. One tool 
which can be used to 'join' polygonal blocks together is the Hadamard product. (For 
a full description, see [16].) This is particularly useful when a few blocks need to be 
joined in a non-recurring manner. We use the Hadamard product in calculating most 
of the generating functions in this paper. 

Consider two series, f(t) = ^2 n f n t n an d git) = ^2, n g n t n , denote by © t the 
Hadamard product with respect to t, then 

f(t)Q t g(t) = J2fn9nt n . (14) 

n 



Exact generating function for 2-convex polygons. 



-1 -4 -4 

x y q 



Figure 4. An example of the action of the Hadamard product. 

The restricted Hadamard product with respect to t is defined as 
f(t) Q t g(t)\t=i = Y,^9n = ^-jf{t)g{l/t)j. 



(15) 



For notational convenience, we will refer to this as a Hadamard join (over t), or simply 
a 'join'. 

Hence, the Hadamard product over a given variable, say s, is the operator which 
'joins' generating functions by matching the perimeters enumerated by s. This is 
equivalent to joining the polygons by matching the edges of the respective polygons, 
such that they overlap, and then removing the overlapping edges, forming a single, 
larger polygon. We can also match the edge column, rather than the edge perimeter. For 
example, if we enumerate staircase and stack polygons according to their right perimeter, 
total perimeter and area, we can join them, making the neighbouring columns overlap. 
Making the transformation s i— > s/yq and dividing by x, so that the overlapping column 
is not double-weighted, we form unimodal polygons, as in figure 4. 

4-2.2. Hadamard arithmetic It is straight-forward to show the following properties of 
restricted Hadamard products [16]. 



It is distributive: 

f{t) Qt (g(t) + h(t)) = f(t) Q t g(t) + f(t) Q t h(t). 



(16) 



It follows the product rule: 
d_ 

ds 



if is, t) Q t g(s, t)) = (—f(s, t)) Q t g(s, t). + f(s, t) Q t (—g(s, t)) (17) 



It evaluates simply at poles: 
f(t) Gt 1 



1 — at 
t k k\ 
(1 - «<)*+! 



/(«) 



t=a 



(18) 
(19) 



We showed in [9] how u and v can be used to simply express the generating functions 
for staircase polygons with fixed steps in the corner. The generating function u a v b counts 



Exact generating function for 2-convex polygons. 



10 




Figure 5. The form of a 1-staircase polygon, which may be formed by inserting an 
indent factor next to a distinguished step. 



staircase polygons that start with a horizontal steps and end with b vertical ones (or 
vice versa). This allows us to enumerate parts of polygons with specific sequences of 
steps along the joins. Moreover, because we can express the generating functions for all 
factors of almost-convex polygons as the quotient of polynomial expressions of u and v, 
evaluating the joins generally becomes straight-forward by separating the poles in the 
denominator using partial fractions. 



4-3. Distinguishing steps to insert indents 

One way of inserting indents in convex polygons is to distinguish a step for the 
location of the indent and then make the appropriate adjustment to the generating 
function. We therefore factorise the polygon at the distinguished step by extending a 
line perpendicular to it. When an indent is joined to a staircase factor (as per figure 5), 
the adjustment required for the indent is independent of its location. 

We distinguish between the walk forming the indent and the rest of that side of 
the polygon. In the case of a vertical indent (as depicted in the figure), the indent 
starts with the vertical step at the same height as the distinguished step preceding 
the non-directed vertical step(s), and includes all steps up to (but not including) the 
distinguished step. We refer to it as the indent factor. We note that a single, m-deep 
indent factor is therefore in the form of a pyramid. 

As mentioned above, the generating function for staircase polygons with a (resp. 
b) fixed horizontal (resp. vertical) steps in one corner is u a v b . As each extra fixed step 
along the factorisation line contributes u to the generating function, the contribution 
of the indent factor to these otherwise staircase polygons (which are called 1-staircase 
polygons) is therefore u 2 /(l — u) 2 . We therefore define the 1-deep indent generating 
function to be 

u 2 v 2 
J =(I^' with J = J ^ = (T^ ( 2 °) 

and, for m > 1, the m-deep indent generating function is 

u 2 v 2 S 2m ~ 2 

l m =X m (x,y)= = , with T m =I m (y,x). (21) 
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(a) The form of 2-staircase polygons (b) The two possible ways of inserting two indents, 
with both indents on their top-left side. 



Figure 6. The form of 2-staircase polygons with both indents on their top-left side. 



The height-independent insertion of indents can be extended to unimodal polygons. 
Whereas staircase polygons can be separated into two halves - one with only positive 
steps, the other with only negative steps - unimodal polygons are defined by the fact 
that their positive horizontal steps occur before the negative ones, and similarly for the 
vertical steps. This means that if the unimodal polygon intersects itself, a staircase 
factor is formed. And so, if we factor a unimodal polygon along the base of a vertical 
indent on the left arc (as per the staircase factorisation shown in figure 5), the indent 
must be a part of a staircase factor to the bottom-left and is therefore enumerated by 
1 m . This leads us to the following proposition. 

Proposition 4.1. The generating function for bimodal m-staircase (resp. m-unimodal) 
polygons that are rooted in the bottom left corner and whose single m-deep indent is 
vertical and on the left arc is Z m y 2 -j^(Q/y), where Q = S/y (resp. xZ). 

This argument can be extended to enumerate unimodal polygons with two distinct 
indents on the same side. To construct such polygons, one may insert both indents at 
a distinguished height, and then mark a second height where we would like the second 
indent. We can therefore try to form the desired polygons by translating the closest of 
the two indents to the second of the distinguished heights. In [12], we point out that if 
the second indent lies below the first, this downward translation may cause the polygon 
to intersect. Furthermore, in the other case, an upward translation will mean that the 
bottom arc of the polygon will not go above the original position of the translated indent. 
We show, however, that the missing polygons in one case are equal in number to the 
extra polygons in the other case. This intriguing fact leads to the following proposition. 



Proposition 4.2. The generating function for 2-staircase (resp. 2-unimodal) polygons 
that are rooted in the bottom left corner and with distinct vertical indents on the left arc 
is 

d d 

y 3 —(l 2 — (Q/y))/2, where Q = S/y (resp. xZ). (22) 
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4-4- Folding Walks and Wrapping Polygons 



4-4- 1- Folding Walks An important notion in deriving exact generating functions for 
almost-convex polygons is that of the so-called 'folding' of DWs. This notion is simply 
a way of describing the combinatorial objects enumerated by the enumerative methods 
developed in [5] . These methods use the half-perimeter operator defined in Section 2 to 
enumerate intersecting convex polygons, and then remove those that intersect. We are 
therefore interested in the enumeration of these intersecting polygons. 

We begin by considering DWs, which are enumerated by 1/(1 — x — y), where x 
(resp. y) counts the horizontal (resp. vertical steps). If we then only consider those with 
an even number of horizontal steps, we can find either a vertex or a series of vertical 
steps that have half of these horizontal steps to the left, and half to the right. We can 
then fold the walk vertically at this point by reflecting all the steps to the left over to the 
right. (The axis of reflection is the vertical line that goes through the half-way point.) 
This walk is now half as wide as it was, and its generating function is therefore 

1 (23) 



E x 



E, r 



(24) 



1 — x — y 

By removing those that have a horizontal step after the fold, we obtain self-avoiding 
walks, which is an example of standard inclusion-exclusion techniques, 

1 — x 
1 — x — y 

Adding a width-one column to its side then gives us stack polygons (that is, sideways 
pyramids) . 

Now, folding vertically as well allows us to make the walk end at its origin, forming 
a polygon. Forcing the polygon to start with a horizontal step and end with a vertical 
one (or vice versa) then gives the following generating function for (possibly intersecting) 
unimodal polygons: 



E 



xy(l - x)(l - y) 



(25) 



1-x-y 

Finally, by distinguishing a horizontal step after which we fold (rather than choosing 
the half-way point) forces a second fold on the other side of the walk in order to make 
its width equal to half its original width. This means that the resulting polygon does 
not necessarily visit any corner of the MBR, and we are left with (possibly intersecting) 
convex polygons. The resulting generating function is 



E 



xy(l 



x) 



yf 



• (26) 

which is the d = 2 case of Bousquet-Melou and Guttmann's multi-dimensional result 
(Lemma 2.2, [5]). 

Folding walks is therefore a simple way of enumerating intersecting polygons (or 
factors of polygons, for that matter) with convexity requirements. This will allow us, in 
the following section, to enumerate 1-convex polygons in a direct, closed-form expression 
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n 



(a) The form of the folded top (b) The form of the wrapped polygon, 

factor. 

Figure 7. The action of wrapping when the top factor is narrower than the length of 
the join, d. 

that is combinatorially interpretable. This provides us with an example of the methods 
required for the enumeration of the various sub-classes of 2-convex polygons. However, 
first we need to be able to fold one factor within a polygon (that is joined directly to 
another factor) without our methods breaking down. And for this we need "wrapping"... 



Jf..Jf..2. Wrapping Polygons 'Wrapping' refers to folding a single factor of a polygon 
that has been constructed by joining multiple factors together. As an example consider 
1-unimodal polygons formed by joining a staircase bottom factor to a unimodal top 
factor and an indent factor. The total height of the polygon is given by the sum of the 
heights of the top and bottom factors. The total width of the polygon is measured by 
adding the width of the top factor to the width of the bottom factor that lies to the left 
of the top factor. In the 1-unimodal case, the generating function for the bottom factor 
is therefore Xv(u/x) d , where d is the length of the join. 

In this case, the unimodal factor is enumerated as a folded walk, as per the previous 
section. However, we will sometimes fold the fixed steps of the polygon, leaving a chain 
of double-bonds of fixed steps around the fold, as shown in figure 7(a). If the join is of 
length d, then there are d fixed horizontal steps in each of the top and bottom factors 
that are identified, but are then removed and do not form part of the polygon. And 
so, if there are fewer than d horizontal steps in the rest of the top unimodal factor, 
as depicted in the figure, then the contribution to the polygon is a pyramid of width 
2n — d, with a weight of x n . Importantly, although those fixed steps are not part of 
the polygon itself, they do contribute to the weight. This is because the width of the 
polygon here is given by the top factor. When the E operator folds the top factor, in 
order to keep the bottom factor joined to the top factor, it must therefore also be folded. 
We say that we have wrapped the bottom factor. We can therefore see that the fixed 
steps along the join that were folded to form double-bonds and whose width is counted 
in the resulting generating function is the projection of the wrapped part of the bottom 
factor. This is shown in figure 7(b). We therefore need not make any adjustment, as the 
width is enumerated correctly. Finally, the required '1-unimodal' polygons are obtained 
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(a) The inclusion case. (b) When the indent extends to the left. 

Figure 8. Generalising Lin's factorisation to enumerate 1-convex polygons. 



by translating the fixed vertical step below the join to the right edge of the polygon (as 
shown in red in figure 7(b)) to ensure that the polygons are self- avoiding. 

In conclusion, without needing to make any extra adjustment, wrapping allows us 
to enumerate almost-unimodal polygons with a single indent in the left side by only 
enumerating staircase polygons joined to unimodal ones. This is because the wrapping 
action also creates polygons that are composed of a unimodal bottom factor joined to 
a pyramid top factor. 



5. The 1-convex generating function 



As an example of the above wrapping technique, we derive the generating function of 1- 
convex polygons with their indent on the top arc. It provides a much simpler derivation 
of the result than the method used in [9]. 

We begin by adopting Lin's factorisation of 1-convex polygons (see [7]) by extending 
a line along the base of the indent, as in figure 8. We enumerate these polygons by 
following the inclusion-exclusion argument of [5], enumerating all the required polygons, 
including those that intersect, and then excluding those that intersect. The factorisation 
gives top and bottom unimodal factors that can be enumerated as walks that may be 
wrapped such that the bottom factor extends furthest to the right. 

The generating function of polygons of the form shown in part (a) of the figure can 
be expressed as 

1 / , . " ./ //- / y 
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The length of the join is taken to be n. The term in the first line is the generating 
function for the unimodal top factor, with a base of at least length n. The term 
x n+1 y 2 /{l — x — y) enumerates DWs that are folded by the E operator to form the 



(27) 
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unimodal factor. We note that due to the vertical symmetry, when this top factor is of 
height one, some polygons may be double-counted, giving the term (1 +y/(l — x)), the 
y/(l — x) forcing the top factor in the symmetric case to be of at least height two. The 
second line enumerates the bottom factor and the indent. Again, using an inclusion- 
exclusion approach, we fold a walk to form the unimodal polygon and then exclude the 
intersecting cases. As the width along the join is already counted in the top factor, we 
adjust for the length of the join with the term (x/ (1 — y)) 71 ' 1 - We then fix the horizontal 
steps of the indent factor, which cannot be folded, giving the term 

( x 2 /(i-y) V 

\l-x 2 /{l-y)) ' 

rather than the expected x 2 /(l — x) 2 , which can be folded. We recognise the last term 
in both of the first two lines as the exclusion cases when the polygons intersects in the 
top-right or bottom-left corner. The last term comes from the possibility of intersection 
in the bottom-right corner. 

From the wrapping principle outlined in the previous section, when the top factor 
has the form of a pyramid, the fixed steps along the join may have been folded. This 
then wraps the bottom factor, making it convex in shape, such that it extends further 
to the right than the top factor. 

We finish by expanding the brackets, so that we may evaluate the sums and write 
the expression in a closed form. However, the indent may extend further to the left 
than the bottom factor, as shown in part (b) of the figure, and therefore requires an 
adjustment term. And so, moving all terms incorporating n into the E operators and 
expanding the brackets in the summand, we can then complete the summation. This 
then allows us to add the term x*/(l — x*) that enumerates the indent for the required 
adjustment to give the following expression for the generating function: 
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(28) 



6. Derivation of the 2-convex generating function 

We generalise Lin's factorisation of 1-convex polygons by extending a factorisation line 
along the base of each indent. This allows each case to be enumerated by joining factors 
along these lines. When the indents are in the same direction, there are three main 
factors plus two indents. When the indents are in different directions, we divide the 
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lattice into four quadrants. We say the quadrant in the top-right is the first, and order 
the remaining quadrants in an anti-clockwise fashion. 

To break-up the problem into enumerable parts, similar to Lin, we adopt a 'divide 
and conquer' approach and classify sub-classes of 2-convex polygons according to the 
relative direction and position of the indents (that is, which side they lie on). We obtain 
the generating functions of symmetric classes by reflection and rotation. Without loss of 
generality we assume that one indent is vertical and on the top-left side. This leaves us 
with one of nine cases: Firstly the indent can have depth two, secondly when there are 
two indents each of depth one, we have to consider eight combinations of the direction 
and location of the second indent. However, the two cases where the second indent is 
in a different direction and on one of the adjacent sides are equivalent after a rotation. 
We are therefore left with eight distinct cases to evaluate. 

There is a possible ambiguity when an indent factor is adjacent to the MBR. For 
example, in the case where there is a vertical indent on the base of the polygon as well 
as on the top-left side, the indent on the base may be considered as on either side of 
the bottom arc. We arbitrarily chose that such cases be enumerated by the class whose 
indented sides are closest together. This example is therefore enumerated by the case 
where the indents are on adjacent sides, and not on opposite sides of the polygon. We 
now briefly consider these eight cases in turn. 

6.1. Case 1: a single 2-deep indent 

Almost-convex polygons with a single indent we refer to as 'bimodal' due to the 
two modes (in the same sense as 'unimodal') adjacent to the indent. In section 5 
we enumerated 1-convex polygons, which are the simplest case of bimodal polygons. 
Bimodal 2-convex polygons, which have a single 2-deep indent, can be enumerated 
mutatis mutandis. For the generating function, see Section 2 of [12]. 

6.2. Case 2: indents in the same direction on the same side 

When the indents are on the same side (on the top- left), we join unimodal top and 
bottom factors to a staircase factor in the middle. This implies that the top factor 
extends furthest to the right. (See figure 9.) This creates three blocks separated by two 
factorisation lines. Each pair of blocks are joined by matching the top and base edges 
along the factorisation line dividing them. Since there is an indent factor in the top 
block, horizontal steps must be added on either side of the indent so their perimeters 
match. In this way the length of the join, counted by the parameter s, is equal to 
the top perimeter of the middle staircase factor. Whenever the join is defined in this 
way, in order to enumerate the indent, a term s 2 /(l — s) 3 is included on one side of 
the Hadamard product in the expression for the generating function. Indeed, this term 
appears in the majority of expressions involving 2-convex polygons. Fortunately, such 
joins can be re-expressed in terms of the first three moments of the generating function 
with the indent omitted. 
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Figure 9. The form of a 2-convex polygon with both indents on the same side. 



In order to match lengths correctly along the bottom join (defined as the length of 
the base of the middle staircase factor and enumerated by the parameter t) we append 
horizontal steps to the bottom factor. In this example, the length of both joins are 
defined as the length of the middle factor along the factorisation lines, as indicated by 
the shaded regions in figure 9. The choice for the definition of the join is determined 
by the generating functions of the top and bottom factors, which are simple rational 
expressions in terms of s and t. This is generally simpler than trying to adjust the 
middle factor to match the other factors. 

In using standard inclusion-exclusion techniques for enumerating the top and 
bottom unimodal factors, wrapping ensures that the cases where the middle or bottom 
factor extends furthest to the right are included. This means that the bottom factor may 
intersect (forming a unimodal loop) to the bottom-right. This is a powerful example of 
the robustness of the wrapping methodology - we fold the top factor, wrapping not just 
the middle factor, but the bottom factor as well, such that it can be the one to extend 
furthest to the right. 

We now complete the enumeration to serve as an example for the following cases. 
We only give some details here, as there are dozens of very complicated formulae in all, 
each using the same principles in their derivation. We outline each of the remaining 
cases to the minimum extent that would be necessary to reproduce the results. We 
break the problem into two parts: when the indents are at the same height and when 
they are not. The latter case is enumerated by the following expression, and the former 
follows mutatis mutandis. 
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where S is the generating function of staircase polygons by base and top perimeter. 
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(a) The left-side indent is (b) The left indent is below the (c) Both indents are on the 
above the corner indent. corner indent. same side. 



Figure 10. The form of 2-convex polygons with two vertical indents on adjacent sides 
of the top arc. 

The expression to the left of the join in the first line enumerates the top factor, 
with the indent enumerated by s 2 / (1 — s) 3 . The first term in the E operator enumerates 
the (possibly intersecting) unimodal folded walks, with the term s/(x — s) counting 
the fixed steps along the join (weighted by xs, but divided by x 2 to adjust for the 
width already enumerated by the middle factor). The term y/(l — x) ensures that it 
is at least of height one. The inclusion-exclusion principle then lends us to exclude 
the intersecting case enumerated in the second term. The second line enumerates the 
bottom factor, obtained by folding a stack polygon (a reflection in the horizontal axis). 
This ensures that the indent is not wrapped. The first term is the part of the stack 
polygon under the join. The last term is the part under the indents. The second term 
enumerates the part of the polygon that extends to the left, including the possibility 
that the bottom factor is a pyramid and that the indent factor extends furthest to 
the left. Finally, the terms in the last line enumerate the exclusion cases. The first 
term counts the polygons that intersect in the bottom-left corner. The 2-unimodal 
factor is enumerated by distinguishing the heights where the indents are placed, thus 
explaining the derivatives. The center term of the derivative is the generating function 
for unimodal polygons of at least height two. The last term counts polygons intersecting 
in the bottom-right corner. This term is obtained by summing the generating functions 
for each possible configuration of the indents, using the known generating function for 
staircase polygons with fixed steps in the corner. 

6.3. Case 3: indents in the same direction on adjacent sides of the same arc 

In this case the top factor must be a pyramid. We break the enumeration into two parts 
depending on whether both indents are at the same height or at different heights. For 
the latter case, illustrated in figure 10, we assume that the top-left indent is higher than 
the top-right one. The symmetrical case can be obtained by reflection. 

In part (c) of the figure, we reproduce the situation from Case 2 (the indents are 
on the same side) where the top factor is a pyramid (such that the middle factor is 
wrapped). Reflecting the indented top pyramid factor vertically we recover the form of 
the polygons shown in part (a). Alternatively, drawing the factorisation line under the 
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(a) The top indent is below the (b) The top indent is level (c) The level indents are 

bottom indent. with the bottom indent. inter-weaved. 




(d) The top indent is above (e) The middle factor is further (f) The indent factor is farthest 
the bottom indent. to the right than the top factor, to the right. 



Figure 11. The form of 2-convex polygons with indents on opposite sides. 



top indent, nipping both the indent and the top factor gives the polygons shown in part 
(b). 

We can therefore use the expression for the generating function for Case 2, with a 
slight modification to the first term, which becomes 
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where P(x,y) is the pyramid generating function. The pyramid term comes from the 
possibility that the indent extends furthest to the right. The last term in the E operator 
excludes those cases from Case 2 that were not wrapped and therefore did not have the 
pyramid top factor now required. 



6-4- Case 4-' indents in the same direction on opposite sides 

We divide this case into three parts according to whether the top-left indent is above, 
next to or below the bottom-right indent (see figure 11.) When the top indent is above 
the bottom one, and either the top or bottom factor is of height one, the indents can be 
considered to be on adjacent sides, and we do not include them in this case. Therefore, 
all of the above calculations must be done for top and bottom factors of height at least 
two. 

When the top indent is below the other (see part (a) of the figure), the enumeration 
does not require wrapping and has no complications; it can be obtained by simply joining 
the three appropriate factors together. When the indents are next to one another (see 
part (b) of the figure), the enumeration is even simpler. However, when the indents are 
level the top indent may be either to the left or the right of the bottom-right indent. 
The latter case means that the indents interweave (see part (c) of the figure). 
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Figure 12. The form of 2-convex polygons with indents in the same direction on 
adjacent sides of opposite arcs. 

Now consider the situation where the top-left indent is above the other one. When 
wrapping, in order to keep the polygon self-avoiding, we usually translate the fixed 
vertical step from below the join to the fold (see figure 7(b)). However, when wrapping 
the bottom factor as well we need to make an adjustment to the length of the bottom 
join (see part (e) of the figure). Moreover, when the middle factor is a pyramid, it is 
possible that the indent extends furthest to the right (see part (f) of the figure). 

This case presents a complication not met earlier in our enumeration. Proceeding 
as usual would involve simultaneously joining an indent and unimodal factor to both the 
top and bottom of a staircase middle factor. When the middle factor extends furthest 
either to the left or right wrapping generates these polygons. However, this requires 
wrapping both sides of the middle staircase factor independently, which is not possible 
as the calculations diverge. (For further discussion, we refer to Section 7.2.1 of [12].) 
We can therefore only use wrapping on one side, and must break up the calculation of 
the bottom factor into two parts depending on whether it is unimodal or pyramid. 

Finally, note that when the middle factor extends furthest to both the left and 
right, it is convex in form. This requires the enumeration of convex polygons according 
to both base and top-most horizontal segment. This can be achieved by adopting a 
'divide and conquer' approach and joining pyramid factors to unimodal ones, or by 
solving recurrence relations. (The generating function by perimeter and area, solved by 
the aforesaid recurrence relations, already exists in the literature [4].) 

6.5. Case 5: indents in the same direction on adjacent sides of opposite arcs 

This case is similar to the previous one, when the bottom factor was a pyramid. 
Reflecting the bottom factor and indent the correct form is produced as shown in 
figure 12. We can proceed by joining a unimodal top factor to a unimodal middle 
factor, which in turn is joined to a pyramid bottom factor. Wrapping again generates 
the cases where the bottom factor extends furthest to the right. The enumeration of 
the different parts of this case follows the same approach as the previous case, except 
that the top and bottom factors may be of height one. Also, the interwoven case is 
somewhat different, as the top and bottom factors are joined directly. 
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(a) The indents form a concave region in the (b) The indents form a convex region in the 
corner. corner. 

Figure 13. The form of 2-convex polygons with indents on the same side in different 
directions. 



In defining these eight cases we have assumed that the top indent is on the top-left 
side. When multiplying by two to obtain the generating function for the symmetric case 
(where the top indent is on the top-right side), we double-count polygons whose top and 
bottom indents are adjacent to the MBR (i.e. they form the topmost segment and the 
base). We must therefore adjust for this case when adding up the generating functions 
at the end. 

6.6. Case 6: indents in different directions on the same side 

When the indents are in different directions and on the same side they form either a 
locally convex or concave region, depending on their order. This is shown in figure 13. 

First consider the locally concave case. Aside from the walk between the two 
indents, the polygon can either enter the first or the third quadrant, but not both. 
If the polygon does not enter the third quadrant then it has only three factors: two 
unimodal ones joined to the staircase factor with the indents. When it passes through 
the third quadrant, the indents and the connecting walk of the first quadrant are joined 
to unimodal factors in the second and fourth quadrants. These are in turn joined by 
a directed walk in the third quadrant. Making the expressions for the second and 
fourth quadrants factor into two parts, one independent of the horizontal join, the other 
independent of the vertical join, allows the calculation of the generating function. This 
is a good example of how wrapping works in both directions to generate all the required 
polygons. 

When the indents form a convex region, the polygon does not enter the first 
quadrant, except in the case where the indents intersect and there are two steps in 
the first quadrant (as shown in the second diagram in part (b) of the figure). The two 
cases are evaluated separately, but the second is a simple version of the first. A unimodal 
factor and an indent in each of the second and fourth quadrants are joined to a staircase 
factor in the third quadrant. The evaluation is straightforward and wrapping gives the 
remaining polygons. 
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(a) The top indent is above (b) The left indent is next (c) The left indent is 

the one in the corner. to the vertical one. below the other. 




(d) The indent factor is adjacent to the vertical (e) The polygon does not enter 

indent. the third quadrant. 



Figure 14. The form of 2-convex polygons with indents in different directions on 
adjacent sides. 

6. 7. Case 7: indents in different directions on adjacent sides 

We enumerate this case by breaking it up into parts classified by the relative height 
of the indents. We say that the indents are next to each other if the holes formed 
overlap in height. If the vertical projection of the humps formed by the indent overlap, 
we say that they are adjacent. And so, the horizontal indent is either 'above', 'next 
to', 'adjacent to' or 'below' the vertical indent (see figure 14). The cases where the 
polygon extends furthest to the left in the third quadrant are all evaluated using the 
usual inclusion-exclusion and wrapping arguments. 

When the horizontal indent is above the vertical one, the polygon may extend 
furthest to the left in the second quadrant (as shown in part (e) of the figure). We need 
to enumerate 1-unimodal polygons whose indent is in the corner according to their base. 
This can be done by joining an almost-pyramid polygon to a unimodal one. These are 
then joined to the bottom factor together with an indent factor, noting that the indent 
may extend further to the right than the bottom factor. 

6.8. Case 8: indents in different directions on opposite sides 

The final case is enumerated in a similar fashion to the previous case. The sub-classes 
defined by the relative heights of the indents are shown in figure 15. The indents need 
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(a) The top indent is above the one in the corner. (b) The left indent is next 

to the vertical one. 




(c) The indent factor is ad- (d) The left indent is (e) The polygon does not 

jacent to the vertical indent. below the other. enter the third quadrant. 

Figure 15. The form of 2-convex polygons with indents in different directions on 
adjacent sides. 

not lie along the MBR and thus the factors in the same quadrant as an indent must be 
of height or width at least two. 

There is, however, one special case to consider, depicted in part (e) of the figure, 
namely when the polygon does not enter the third quadrant. This is an interesting 
case because the horizontal and vertical joins (over s and t respectively) must be done 
simultaneously. We cannot proceed is the usual way and are forced to re-evaluate our 
approach and come up with a new way to evaluate such constructions. 

We solve this problem by generating the factor in the first quadrant as a staircase 
factor and wrapping it along the factorisation lines to form a convex factor. This 
is achieved by adding double-bonds joined to the staircase factor next to the indent 
factors, so that it wraps back along the double bonds, forming a convex polygon. The 
length of the horizontal (vertical) join is counted by s (t). Each horizontal double-bond 
has two steps, contributing s 2 to the weight of the join, as well as a weight of 1/x for 
the fact that it wraps the staircase back, reducing its width by one. This means that 
the expression for the bottom factor in the join includes the term s 2 / (x — s 2 ). 

Other expressions including s 2 rather than just s have been evaluated preciously by 
simplifying the other side of the join and re-expressing it in terms of derivatives. Now 
what is on the other side of the join is the staircase factor, which cannot be simplified as 
it has to be joined to both the bottom and left factors at the same time. We are therefore 
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forced to simplify the terms including s 2 or t 2 . However since s 2 / (x — s 2 ) = E x [s/ (x — s)], 
the expression inside the E operator can be changed into a form we can evaluate simply. 

The remaining cases are obtained using the above methodology by wrapping a 
directed walk in the first quadrant. When the top indent is above the bottom one the 
unimodal factor in the third quadrant must be of width and height at least two . 

6.9. The 2-convex generating function 

In the above derivations, the direction and position of one of the indents was chosen 
arbitrarily, such that the direction and position of the other determined which sub-class 
the polygon belongs to. Therefore, when the two indents are in different directions, 
the cardinality of the set of possible combinations of direction and location for the 
fixed indent is four. When the indents are in the same direction, both directions 
are enumerated by the generating functions. If the indents are on the same edge, 
the cardinality is two. If they are opposite, it is one. We therefore obtain the 
generating function for 2-convex polygons by summing the results, multiplying each 
term by the cardinality of the class it counts. This gives the following generating 
function for 2-convex polygons (available in a Maple/Mathematica friendly format at 
http:/ /www. ms.unimelb.edu.au/~iwan/polygons/series/2-convex-anisotropic.txt): 
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3925x 2 - 2349x 3 + 873x 4 - 637x 5 -411x 6 + 504x 7 + 103x 8 + 10x 9 )y 10 + (-330+ 1422x - 1797x 2 + 
434x 3 + 47x 4 + 484x 5 - 100x 6 - 160x 7 - 22x 8 )y n + (165 - 603x + 632x 2 - 4x 3 - 180x 4 - 114x 5 + 
88x 6 + 18x 7 )y 12 - (1 - x) 2 (55 - 67x - 18x 2 + 27x 3 + 7x 4 )y 13 - (1 - x) 3 (-ll + x 2 )y 14 - (1 - x) 3 y 15 

and B ( 2 c) = 4(l-x) 26 x 4 -4(l-x) 24 x 4 (30 - 20x+7x 2 )y+4(l-x) 22 x 2 (4-12x+433x 2 - 580x 3 + 
368x 4 -128x 5 +21x 6 )y 2 -4(l-x) 20 x 2 (100-360x+4214x 2 -8192x 3 + 7616x 4 -4297x 5 + 1528x 6 - 
350x 7 +35x 8 )y 3 +4(l-x) 18 (l-12x+1252x 2 -5226x 3 +32426x 4 -76437x 5 +93156x 6 -70591x 7 + 
36061x 8 - 12917x 9 + 3217x 10 -532X 11 + 35x 12 )y 4 - (1 -x) 16 (104- 1296x+42936x 2 - 198752x 3 + 
859020x 4 -2164230X 5 +3180634x 6 -3007804x 7 +1964426x 8 -927123x 9 +322834x 10 -83390x 41 + 
15842x 12 - 1961x 13 + 84x 14 )y 5 + (1 -x) 14 (1300- 16868x + 289112x 2 - 1412308x 3 + 5065708x 4 - 
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12664938x 5 +20907814x 6 -23411754x 7 + 18540432x 8 -10767204x 9 +4716295x 10 -1577652x n + 
403402x 12 -79358x 13 + 12073x 14 - 1114x 15 +28x 16 )y 6 - (1 -x) 12 (10400- 140832x+ 1632176x 2 - 
8077604x 3 + 26714752x 4 - 64771716x 5 + 113604188x 6 - 143655868x 7 + 132836186x 8 - 
91595522x 9 +48102042x 10 - 19641626X 11 +6301753x 12 - 1574826x 13 +304364x 14 -47486x 15 + 
6143x 16 - 400x 17 + 4x 18 )y 7 - (1 - x) 10 (-59800 + 846768x - 7892920x 2 + 38780872x 3 - 
125239164x 4 + 297723354x 5 - 537008734x 6 + 736108516x 7 - 767512896x 8 + 611776002x 9 - 
375779397x 10 + 180183183X 11 - 68757620x 12 + 21225217x 13 - 5236463x 14 + 983039x 15 - 
136742x 16 + 16579x 17 - 1952x 18 + 78x 19 )y 8 + (1 - x) 8 (-263120 + 3902272x - 32692880x 2 + 
159495728x 3 - 518222904x 4 + 1235851899x 5 - 2273834564x 6 + 3286427610x 7 - 3744105010x 8 + 
3360169186x 9 -2372247082x 10 + 1316981454x n -578123374x 12 +204616292x 13 -60356681x 14 + 
15025753x 15 - 2968864x 16 + 408498x 17 - 36191x 18 + 3653x 19 - 450x 20 + 7x 21 )y 9 + (1 - 
x) 6 (920920 - 14323848x + 115445704x 2 - 565821908x 3 + 1881114592x 4 - 4595055622x 5 + 
8689038821x 6 - 13096803662x 7 + 15935786206x 8 - 15697620474x 9 + 12485332431x 10 - 
7965917640X 11 + 4040636749x 12 - 1619517767x 13 + 518745262x 14 - 140030433x 15 + 
34601374x 16 - 7726101x 17 + 1290999x 18 - 115181x 19 + 1171x 20 - 68x 21 + 43x 22 )y 10 + (1 - 
x) 4 (-2631200 + 42965472x - 346781248x 2 + 1732729172x 3 - 5960460908x 4 + 15133786878x 5 - 
29768776966x 6 + 46871962702x 7 - 60257832295x 8 + 63871376160x 9 - 55950160750x 10 + 
40346497476X 11 - 23700589523x 12 + 11139364197x 13 - 4087427426x 14 + 1145559467x 15 - 
251575315x 16 + 52651169x 17 - 13478667x 18 + 3310337x 19 - 518486x 20 + 32936x 21 + 909x 22 + 
2x 23 + 3x 24 )y n - (1 - x) 3 (-6249100 + 100966228x - 786008036x 2 + 3790621220x 3 - 
12634194004x 4 + 31120468774x 5 - 59333283779x 6 + 90481426644x 7 - 112836337344x 8 + 
116619494723x 9 - 100526269572x 10 + 72257250633X 11 - 42935196678x 12 + 20646926051x 13 - 
7717297655x 14 + 2081066620x 15 - 344551302x 16 + 19708982x 17 - 190245x 18 + 2604586x 19 - 
108 9 9 25x 20 + 167 9 90x 21 - 53 25x 22 - 8 19x 23 + 5x 24 )y 12 - (1 - x) 2 (12498200 - 200432672x + 
1523814896x 2 - 7173568104x 3 + 23394202064x 4 - 56432397385x 5 + 105258925336x 6 - 
156756419659x 7 + 190737577732x 8 - 192647570536x 9 + 163201860685x 10 - 116641434399X 11 + 
70265808073x 12 - 35189781755x 13 + 14129858643x 14 - 4206300270x 15 + 752961328x 16 + 
2762280x 17 -41077210x 18 +8781488x 19 + 172960x 20 -372435x 21 +59206x 22 -2116x 23 -89x 24 + 
llx 25 )y 13 -(l-x)(-21246940+339297288x-2543202516x 2 +11796257160x 3 -37956243060x 4 + 
90389633588x 5 - 166278246916x 6 + 243693863948x 7 - 291099521233x 8 + 288163147560x 9 - 
239402103854x 10 + 168782468152X 11 - 101938448323x 12 + 52836781903x 13 - 23089813175x 14 + 
8076150924x 15 - 2006509792x 16 + 229306444x 17 + 52941397x 18 - 29443771x 19 + 5233574x 20 - 
137 3 08x 21 - 96 84x 22 + 130 23x 23 - 102x 24 + 33x 25 )y 14 + (-30904640 + 493047872x - 
3671274848x 2 + 16909765768x 3 - 54076869168x 4 + 128034624944x 5 - 233957204524x 6 + 
339860131708x 7 -401161475162x 8 +390960587591x 9 -318456383150x 10 + 219480309921X 11 - 
130108934416x 12 + 67655799311x 13 - 31165964099x 14 + 12455845878x 15 - 4042728520x 16 + 
929084638x 17 - 93195098x 18 - 24059794x 19 + 11705293x 20 - 1983652x 21 + 86582x 22 + 
23 5 47x 23 - 3 151x 24 - 68x 25 + 5x 26 )y 15 + (38630800 - 579105408x + 4018988224x 2 - 
17153382704x 3 + 50579891852x 4 - 109880934764x 5 + 183275392200x 6 - 241692656260x 7 + 
257400864147x 8 - 224539648706x 9 + 161936707748x 10 - 97616424050X 11 + 50386570178x 12 - 
23295420511x 13 + 10082660716x 14 - 4042728520x 15 + 1381926138x 16 - 352763968x 17 + 
51556603x 18 + 1631397x 19 - 2559687x 20 + 547928x 21 - 54439x 22 - 611x 23 + 432x 24 - 31x 25 )y 16 + 
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(-41602400+585761792x-3791089760x 2 + 15002579968x 3 -40793586584x 4 +81262920502x 5 - 
123537868000x 6 + 147559428776a; 7 - 141331365760x 8 + 109736530568x 9 - 69193200719x 10 + 
35421383132X 11 - 15013860403x 12 + 5709460646x 13 - 2235816236x 14 + 929084638x 15 - 
352763968x 16 + 103409188x 17 - 19757856x 18 + 1570918x 19 + 189358x 20 - 71908x 21 + 
11883x 22 - 268x 23 + 77x 24 )y 17 + (38630800 - 510758608x + 3084598320x 2 - 11324232712x 3 + 
28395037220x 4 - 51826962892x 5 + 71682098798x 6 - 77353955608x 7 + 66429125667x 8 - 
45706907743x 9 + 24906199156x 10 - 10403200527X 11 + 3174037717x 12 - 706359558x 13 + 
176365047x 14 -93195098x 15 +51556603x 16 -19757856x 17 +4723092x 18 -605060x 19 +49569x 20 - 
3395x 21 - 1253x 22 - 143x 23 )y 18 + (-30904640 + 383692992x - 2163792656x 2 + 7374341864x 3 - 
17053687520x 4 + 28497776084x 5 - 35790327136x 6 + 34800116028x 7 - 26729224736x 8 + 
16280988977x 9 -7648226004x 10 +2524992746X 11 -406853569x 12 -93698188x 13 +82385168x 14 - 
24059794x 15 + 1631397x 16 + 1570918x 17 - 605060x 18 + 66554x 19 - 4454x 20 + 1819x 21 + 
202x 22 )y 19 + (21246940 - 247772008x + 1305971612x 2 - 4135556872x 3 + 8823936024x 4 - 
13492690928x 5 + 15356870874x 6 - 13416475166x 7 + 9196907145x 8 - 4978222336x 9 + 
2052527732x 10 - 556811827X 11 + 29183400x 12 + 58467226x 13 - 34677345x 14 + 11705293x 15 - 
2559687x 16 + 189358x 17 + 49569x 18 - 4454x 19 - 1368x 20 - 110x 21 )y 20 + (-12498200 + 
136980272x - 675634696x 2 + 1990453776x 3 - 3921712444x 4 + 5486790816x 5 - 5649513210x 6 + 
4421089360x 7 - 2695886222x 8 + 1305174655x 9 - 491084981x 10 + 128534739X 11 - 11441768x 12 - 
8063133x 13 +5370882x 14 -1983652x 15 +547928x 16 -71908x 17 -3395x 18 + 1819x 19 -110x 20 )y 21 + 
(6249100 - 64433028x + 297851620x 2 - 817733112x 3 + 1489860872x 4 - 1908408554x 5 + 
1775572130x 6 - 1240442976x 7 + 667928313x 8 - 289317379x 9 + 102644601x 10 - 29961766X 11 + 
6383656x 12 - 977036x 13 - 41224x 14 +86582x 15 - 54439x 16 + 11883x 17 - 1253x 18 + 202x 19 )y 22 + 
(-2631200 + 25559072x - 110940720x 2 + 284446548x 3 - 480234432x 4 + 564112724x 5 - 
474248936x 6 + 294904852x 7 - 138225516x 8 + 5252241 lx 9 - 17224508x 10 + 5578263X 11 - 
1609051x 12 + 492963x 13 - 109107x 14 + 23547x 15 - 611x 16 - 268x 17 - 143x 18 )y 23 + (920920 - 
8444128x + 34492144x 2 - 82799336x 3 + 129893340x 4 - 140307634x 5 + 106798702x 6 - 
59095890x 7 + 23639299x 8 - 7486656x 9 + 2001655x 10 - 644781X 11 + 181503x 12 - 63349x 13 + 
13125x 14 - 3151x 15 + 432x 16 + 77x 17 )y 24 + (-263120 + 2283072x - 8798288x 2 + 19829008x 3 - 
28995448x 4 + 28906719x 5 - 20012600x 6 + 9913898x 7 - 3347134x 8 + 826000x 9 - 124087x 10 + 
29300X 11 - 2853x 12 + 1927x 13 - 135x 14 - 68x 15 - 31x 16 )y 25 + (59800 - 492568x + 1795744x 2 - 
3809804x 3 + 5207920x 4 - 4809154x 5 + 3046581x 6 - 1369614x 7 + 399732x 8 - 78407x 9 + 2224x 10 + 
919X 11 - 834x 12 + lllx 13 + 33x 14 + 5x 15 )y 26 + (-10400 + 81632x - 282336x 2 + 565028x 3 - 
7231 16x 4 + 619222x 5 - 359946x 6 + 148482x 7 - 39609x 8 + 7449x 9 - 326x 10 - lOx 11 + 5x 12 - 
llx 13 )y 27 + (1300 - 9768x + 32140x 2 - 60712x 3 + 72592x 4 - 57298x 5 + 302 1 7x 6 - 11207x 7 + 
2732x 8 - 506x 9 + 43x 10 + 3x n )y 28 + (-104 + 752x - 2360x 2 + 4200x 3 - 4648x 4 + 3305x 5 - 
1506x 6 + 448x 7 - 78x 8 + 7x 9 )y 29 + 4(1 - x) V°. 

7. Summary and outlook 

Following Lin's approach to enumerating 1-convex polygons, one can factorise almost- 
convex polygons by extending lines along the base of all indents. Then, using a 'divide 
and conquer' approach, it is possible to then enumerate the various sub-classes. 
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By looking at the form of the various factors, it is possible to guess what the form of 
the resulting generating functions will be. One can then obtain the generating functions 
by directly enumerating the series to a sufficiently high order and then solving the set 
of linear equations corresponding to the presumed form. 

We presented some techniques that are invaluable in enumerating the factors 
exactly. The most important technique is that of 'wrapping', which allows the generation 
of quite complex objects out of simply enumerable components. It is easily implemented 
when using an inclusion-exclusion approach to enumerating unimodal factors, which 
made it essential in enumerating 2-convex polygons. These techniques have allowed us 
to reduce the enumeration of 1-convex polygons, as well as many sub-classes of 2-convex 
polygons, to a single, combinatorially interpretable expression. 

Going forward, it is not realistic to factorise almost-convex polygons for high 
concavity indices as we have done here, as there will be an exponential growth in the 
number of cases to evaluate. It would be more sensible to restrict the size of each 
indentation first, and then generalise these cases. Eventually, defining operators that 
can add more and more complex indentations in the side of convex polygons and looking 
at the effect on the asymptotic growth of their number seems to be the most appropriate 
path to understand how convex polygons become general SAPs. 
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